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Two different tune measurement systems have been installed in the GSI heavy-ion syn- 
chrotron SIS-18. Tune spectra are obtained with high accuracy using these fast and sensitive 
systems . Besides the machine tune, the spectra contain information about the intensity de- 
pendent coherent tunc shift and the incoherent space charge tune shift. The space charge 
tune shift is derived from a fit of the observed shifted positions of the synchrotron satellites 
to an analytic expression for the head-tail eigenmodes with space charge. Furthermore, the 
chromaticity is extracted from the measured head-tail mode structure. The results of the 
measurements provide experimental evidence of the importance of space charge effects and 
head-tail modes for the interpretation of transverse beam signals at high intensity. 
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I. INTRODUCTION 

Accurate measurements of the machine tune and of the chromaticity are of importance for 
the operation of fast ramping, high intensity ion synchrotrons. In such machines the tune spread 
fiQx,y at injection energy due to space charge and chromaticity can reach values as large as 0.5. 
In order to limit the incoherent particle tunes to the resonance free region the machine tune has 
be controlled with a precision better than AQ ~ 10~ 3 . In the GSI heavy-ion synchrotron SIS- 18 
there are currently two betatron tune measurement systems installed. The frequency resolution 
requirements of the systems during acceleration are specified as 10 -3 , but they provide much higher 
resolution (10 ) on injection and extraction plateaus. The Tune, Orbit and Position measurement 
system (TOPOS) is primarily a digital position measurement system which calculates the tune 
from the measured position [I]. The Baseband Q measurement system (BBQ) conceived at CERN 
performs a tune measurement based on the concept of diode based bunch envelope detection [2]. 
The BBQ system provides a higher measurement sensitivity than the TOPOS system. Passive tune 
measurements require high sensitivity (Schottky) pick-ups, low noise electronics and long averaging 
time to achieve reasonable signal-to-noise ratio. For fast tune measurements with the TOPOS and 
BBQ systems, which both use standard pick-ups, the beam has to be excited externally in order 
to measure the transverse beam signals. 

For low intensities the theory of transverse signals from bunched beams and the tune measure- 
ment principles from Schottky or externally excited signals are well known [3 H] . In intense, low 
energy bunches the transverse signals and the tune spectra can be modified significantly by the 
transverse space charge force and by ring impedances. Previously, the effect of space charge on 
head-tail modes had been the subject of several analytical and simulation studies [SHS]- 

Recently, the modification of transverse signals from high intensity bunches was observed in the 
SIS-18 for periodically excited beams [TU] and for initially kicked bunches [11], where the modified 
spectra was explained in terms of the space charge induced head-tail mode shifts. 

This contribution aims to complement the previous studies and extract the relevant intensity 
parameters from tune spectra measurements using the TOPOS and BBQ tune measurement sys- 



tems. Section II presents the frequency content of transverse bunched beam signals briefly. Sec 



tion III presents the space charge and image current effects on tune measurements and respective 



theoretical models. Section IV report on the experimental conditions and compares the charac- 



teristics of two installations as well as the various excitation methods. Section V presents the 



experimental results in comparison with the theoretical estimates of various high intensity effects. 
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II. TRANSVERSE BUNCH SIGNALS 

Theoretical and experimental work related to transverse Schottky signals and beam transfer 
functions (BTFs) for bunched beam at low intensities can be found in the existing literature 
[HIES EI]- The transverse signal of a beam is generated by the beam's dipole moment 

N 

3=1 

where Xj(t) is the horizontal/vertical offset and Ij(t) the current of the j th particle at the position 
of a pick-up (PU) in the ring. The sum extends over all N particles in the detector. The Schottky 
noise power spectrum as a function of the frequency / is defined as S(f) = \d(f)\ 2 where d(f) is 
the Fourier transformed PU signal. For a beam excited by an external amplitude spectrum G(f), 
the transverse response function is defined as r(f) = d(f)/G(f) [15J. 

If the transverse signal from a low intensity bunch is sampled with the revolution period T s , 
then the positive frequency spectrum consists of one set of equidistant lines 

Q k = Q + AQ k , 2 

usually defined as baseband tune spectrum, where Qq is the fractional part of the machine tune, 
AQk = ±kQ s are the synchrotron satellites and Q s is the synchrotron tune. 

For a single particle performing betatron and synchrotron oscillations, the relative amplitudes 
of the satellites are [3] 

| d k |~| J fc (x/2) I 3 

where x = ^C^m/vo is the chromatic phase, £ is the chromaticity, cp m is the longitudinal oscillation 
amplitude of the particle and rjo the frequency slip factor. J k are the Bessel functions of order k. 
In bunched beams, the relative height and width of the lines depends on the bunch distribution. 
The relative height is also affected by the characteristics of the external noise excitation or by the 
initial transverse perturbation applied to the bunch. In the absence of transverse nonlinear field 
components the width of each satellite with k ^ is determined by the synchrotron tune spread 
5Qk ~ |&|Q s <^j/16. An example tune spectrum obtained from a simulation code for a Gaussian 
bunch distribution is shown in |Fig. 1 . 
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FIG. 1. Baseband tune spectrum from a simulation code for a low intensity bunch. 

III. TUNE SPECTRUM FOR HIGH INTENSITY 

At high beam intensities, the transverse space charge force together with the coherent force 
caused by the beam pipe impedance will affect the motion of the beam particles and also the tune 
spectrum. The space charge force induces an incoherent tune shift Qo — AQ SC for a symmetric 
beam profile of homogeneous density where 

AQsc = ^p R 4 

is the tune shift, I p the bunch peak current, q the particle charge and Eq = "Jquic 2 the total energy. 
The relativistic parameters are 70 and /3q, the ring radius is R and the emittance of the rms 
equivalent K-V distribution is e x . In the case of an elliptic transverse cross-section the emittance 
e x in 



Eq. 4 should be replaced by 



For the vertical plane the procedure is the same, with x replaced by y. The image currents and 
image charges induced in the beam pipe, assumed here to be perfectly conducting, cause a purely 
imaginary horizontal impedance 

Zx = "StKAjToM 2 6 

and real coherent tune shift 

. qIR 2 Z x 



2Q x oPoEq 

For a round beam profile with radius a and pipe radius b the coherent tune shift is smaller by 

2 

AQ C = AQsc than the space charge tune shift. Therefore the contribution of the pipe is especially 
important for thick beams (a ~ b) at low or medium energies. 
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In the presence of incoherent space charge, represented by the tune shift AQ SC , or pipe effects, 
represented by the real coherent tune shift AQ C , the shift of the synchrotron satellites in bunches 
can be reproduced rather well by [B] 

AQ k = - AQsc ± AQc ± V(AQ SC - AQ c ) 2 /4 + {kQ s y 8 

where the sign + is used for k > 0. For k=0 one obtains AQ^q = —AQ C . The above expression 
represents the head-tail eigenmodes for an airbag bunch distribution in a barrier potential [5] with 
the eigenfunctions 

= cos(kTt(j)/(f>b) exp(— ix4>/4>b) 9 

where x is the local transverse bunch offset, x = £,<frb/vo is the chromatic phase, 4>b is the full 
bunch length and rjo is slip factor. The head-tail mode frequencies obtained from |Eq. 8| are shown 
in Fig. 2 . In Ref. [5] the analytic solution for the eigenvalues Eq. 8 is obtained from a simplified 
approach, where the transverse space charge force is assumed to be constant for all particles. This 
assumption is correct if there are only dipolar oscillations. In Ref. [6] it is has been pointed out, 
that in the presence of space charge there is an additional envelope oscillation amplitude. For the 
negative-A: eigenmodes the envelope contribution dominates and therefore those modes disappear 



from the tune spectrum. In Ref. [6] |Eq. 8 has been successfully compared to Schottky spectra 



obtained from 3D self-consistent simulations for realistic bunch distributions in rf buckets. Analytic 
and numerical solutions for Gaussian and other bunch distribution valid for q sc 3> 1 were presented 
in [8]. 

In an rf bucket the synchrotron tune Q s is a function of the synchrotron oscillation amplitude 4>. 
For short bunches Q s corresponds to the small-amplitude synchrotron tune 

2 qVoh\7] \ 
^ s0 ~ 2^m 7 /? 2 c 2 

where Vo is the rf voltage amplitude and h is the rf harmonic number. 

For head-tail modes the space charge parameter is defined as a ratio of the space-charge tune 



shift ( Eq. 4 ) to the small-amplitude synchrotron tune, 



and the coherent intensity parameter as, 



AQ SC 

Qsc = — 11 

<^s0 



Qc = 12 
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FIG. 2. Head-tail mode frequencies as a function of the space charge parameter using |Eq. 8| . The red 
curves represent the result obtained for q c — q sc /10. 
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FIG. 3. Simulation result for g sc = 2. 

An important parameter for head-tail bunch oscillations in long bunches is the effective syn- 
chrotron frequency which will be different from the small-amplitude synchrotron frequency in short 
bunches. For an elliptic bunch distribution (parabolic bunch profile) with the bunch half-length 
4>m = V^&i (rms bunch length 07), one obtains the approximate analytic expression for the longi- 
tudinal dipole tune [17] , 

Using Q s \ instead of Q s q in Eq. [8] shows a much better agreement with the simulation spectra for 
long bunches in rf buckets (see Ref. 

For Gaussian bunches with a bunching factor Bj = 0.3 (Bf = Iq/I p , Iq is the dc current), the 
transverse tune spectra obtained from PATRIC simulations |6j for different space charge factors 
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FIG. 4. Simulation result for q sc = 10. 
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FIG. 5. Head-tail tune shifts and their width obtained from the simulations for q c = 0. The error bars 
indicate the widths of the lines. 



and thin beams (q c = 0) are shown in Figures [TJ [3] and [4j The dotted vertical lines indicate the 
positions of the head-tail tune shifts obtained from |Eq. 8 with Q s = Q s i- For the \ow-k satellites 



there is a good agreement between Eq. 8 and the simulation results. Lines with k > 2 can only 
barely be identified in the simulation spectra. The positions of the satellites for k = 0, 1, 2 together 



with the predicted head-tail tune shifts from Eq. 8 are shown in Fig. 5 . The error bars indicate 
the obtained widths of the peaks in the tune spectra. 

It is important to notice that the simulations for moderate space charge parameters (q sc < 10) 
require a 2.5D self-consistent space charge solver. The theoretical studies rely on the solution of the 
Mohl-Schnauer equation [12J, which assumes a constant space charge tune shift for all transverse 
particle amplitudes. In contrast to the self-consistent results, PATRIC simulation studies using 
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FIG. 6. Head-tail tune shifts and their width obtained from the simulations for q c — 0. 15g sc . The error bars 
indicate the widths of the peaks. 



the Mohl-Schnauer equation gave tune spectra with pronounced, thin satellites also for large k. In 
order to account for the intrinsic damping of head-tail modes [TpOj , which is the main cause of the 
peak widths obtained from the simulations, a self-consistent treatment is required. 

For thick beams (here q c = 0.15q sc , which corresponds to the conditions at injection in the 
SIS-18) the positions of the synchrotron satellites obtained from the simulation are indicated in 
Fig. 6| . A gain, the error bars indicate the widths of the peaks. From the plot we notice an increase 



in the spacing between the k = 0, 1, 2 satellites, relative to the analytic expression. Also the peak 
width for k = 1, 2 does not shrink with increasing q sc . 



The tune spectra obtained for q sc = 3, q sc = 5 and q sc = 10 are shown in Fig. 7 , Fig. 8 



and Fig. 9 . One can observe that for thick beams (here a ~ 0.46) the k = 1 peak remains very 
broad up to q sc = 10. 

This observation is consistent with a very simplified picture for the upper q sc threshold for 
the intrinsic Landau damping of head-tail modes. For a Gaussian bunch profile, the maximum 
incoherent tune shift, including the modulation due to the synchrotron oscillation is 

AQ max = -AQ SC + kQ s 14 

where AQ SC is determined by Eq. |4j The minimum space charge tune shift is (see Refs. [SI [9]) 

AQ m in = -aq sc Q s + kQ s 15 

where a is determined from the average of the space charge tune shift along a synchrotron oscillation 
with the amplitude <j) = 4> m . For a parabolic bunch we obtain a = 0.5. For a Gaussian bunch 
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FIG. 7. Tune spectrum obtained from the simulation for q sc = 3 and q c = 0.15(7, 
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FIG. 8. Tune spectrum obtained from the simulation for q sc = 5 and q c — 0.l5q sc . 

and <f> = 3o";, we obtain a = 0.287. Each band of the incoherent transverse spectrum has a lower 
boundary determined by the maximum tune shift AQmax and an upper boundary determined by 
AQmin- Landau damping, in its very approximate treatment, requires an overlap of the coherent 
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FIG. 9. Tune spectrum obtained from the simulation for q sc = 10 and q c = 0.15g sc . 
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FIG. 10. Head-tail mode frequencies as a function of the space charge parameter. The grey shaded area 
indicates the incoherent band for a Gaussian bunch and k = 2. The red curves represent the head-tail mode 
frequencies obtained for q c — 0.15g sc . 



peak with the incoherent band. The head-tail tune for low q sc can be approximated as 

AQ k = -^(AQ sc + AQ c ) + kQ s 16 

The distance between the coherent peak and the upper boundary of the incoherent band for fixed 
k is 

SQk = Q - aj AQ SC + ^AQ C 17 

For large q sc the head-tail modes with positive k converge towards Qk = —AQ c /2. For a given k 
the mode is still inside the incoherent band if 

k > aq sc - ^q c 18 

holds. In order to illustrate the above analysis, the incoherent band for k = 2 is shown in |Fig. 10 
(shaded area). For q c = the coherent head-tail mode frequency crosses the upper boundary of 
the band at q sc ~ 4.5. For q c = 0.15q sc the k = 2 head-tail mode remains inside the incoherent 
band until q sc ~ 12. For the k = 1 modes the above analysis leads to thresholds of q sc ~ 2 (thin 
beams) and q sc w 6 for q c = 0.15q sc . 



IV. MEASUREMENT SETUP FOR TRANSVERSE BUNCH SIGNALS 

In this section, a brief description of the transverse beam excitation mechanisms as well as 
the two different tune measurement systems, TOPOS and BBQ, in the SIS-18 is given. Further 
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the experimental set-up, typical beam parameters and uncertainty analysis of the measured beam 
parameters is discussed. 



A. Transverse Beam Excitation 



The electronics used for beam excitation consist of a signal generator connected to two 25 W 



amplifiers which feed power to 50f2 terminated stripline exciters as shown in Fig. 11 . Excitation 
types such as band limited noise and frequency sweep are utilized at various power levels to induce 
coherent oscillations. 



1. Band limited noise: 



Band limited noise is a traditionally used beam excitation system for slow extraction in the 
SIS-18. The RF signal is mixed with Direct Digital Synthesis (DDS) generated fractional tune 
frequency, resulting in RF harmonics and their respective tune sidebands. This signal is further 
modulated by a pseudo-random sequence resulting in a finite band around the tune frequency. 
The width of this band is controlled by the frequency of the pseudo-random sequence. Typical 
bandwidth of band limited exciter is ~ 5% of the tune frequency. There are two main advantages 
of this system; first it is an easily tunable excitation source available during the whole acceleration 
ramp and second, the band limited nature of this noise results in an efficient excitation of the beam 
in comparison to white noise excitation. The main drawback is the difficulty in correlating the 
resultant tune spectrum with the excitation signal. 



2. Frequency sweep: 

Frequency sweep (chirp/harmonic excitation) using a network analyzer for BTF measurements 
is an established method primarily for beam stability analysis [4j. However, using this method for 
tune measurements during acceleration is not trivial, and thus the method is not suitable for tune 
measurements during the whole ramp cycle. Nevertheless, this method offers advantages compared 
to the previous excitation method for careful interpretation of tune spectrum in storage mode, 
e.g., injection plateau or extraction flat top. Thus frequency sweep is used during measurements 
at injection plateau to compare and understand the dependence of tune spectra on the type of 
excitation. 
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B. TOPOS 



Following the beam excitation, the signals from each of the 12 shoe-box type BPMs [18] at 
SIS-18 pass through a high dynamic range (90 dB) and broadband (100 MHz) amplifier chain from 
the synchrotron tunnel to the electronics room, where the signals are digitized using fast 14 bit 
ADCs at 125 MSa/s. Bunch-by-bunch position is calculated from these signals using FPGAs in 
real time and displayed in the control room. The spatial resolution is ~0.5 mm in bunch-by-bunch 
mode. Further details can be found in [I]. Hence, TOPOS is a versatile system which provides 
accurate bunch-by-bunch position and longitudinal beam profile. This information is analyzed to 
extract non trivial parameters like betatron tune, synchrotron tune, beam intensity evolution etc. 




FIG. 11. Scheme of the excitation system and the Tune, Orbit and Position (TOPOS) measurement system. 



C. BBQ 



The BBQ system is a fully analog system and its front end is divided into two distinct parts; 
a diode based peak detector and an analog signal processing chain consisting of input differential 
amplifier and a variable gain filter chain of 1 MHz bandwidth. The simple schematic of BBQ 



system configuration at SIS-18 is shown in Fig. 12 and the detailed principle of operation can be 
found in Ref. [2]. 
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FIG. 12. BBQ: Baseband Q measurement system. Diode detectors (a) and signal chain (b). 



D. Comparison of TOPOS and BBQ 



The sensitivity of BBQ has been measured to be ~10-15 dB higher than that of TOPOS under 
the present configuration. The main reason for the difference is the relative bandwidth of the two 
systems and their tune detection principles. In BBQ, the tune signal is obtained using analog 
electronics (diode based peak detectors and differential amplifier) immediately after the BPM 
plates, while in TOPOS position calculation is done after passing the whole bunch signal through 
a wide bandwidth amplifier chain. Even though the bunches are integrated to calculate position in 
TOPOS which serves as a low pass filter, the net signal-to-noise ratio is still below BBQ. Operations 
similar to BBQ could also be performed digitally in TOPOS to obtain higher sensitivity but would 
require higher computation and development costs. TOPOS can provide individual tune spectra 
of any of the four bunches in the machine, while BBQ system provides "averaged" tune spectra of 
all the bunches. Both the systems have been benchmarked against each other. Tune spectra shown 
in Sec. V are mostly from the BBQ system while TOPOS is primarily used for time domain 



analysis, nevertheless this will be pointed out when necessary. 
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TABLE I. Beam parameters during the N 7+ experiment 



Beam/Machine parameters 


Symbols 


Values 


Atomic mass 


A 


238 


Charge state 


q 


73 


Kinetic energy 


Ehin 


11.4 McV/u (measured 


Number of particles 


N p 


1,5, 12 • 10 8 (measured) 


Tune 


Qx 7 Qy 


4.31, 3.27 (set value) 


Chromaticity 


£x i £y 


-0.94, -1.85 (set value) 


Transverse emittance 


e x ,e y (2a) 


45, 22 mm-mrad (measured) 


Slip factor 


V 


0.94 


Bunching factor 


Bf 


0.4 (measured) 


Synchrotron tune 


QsOj Qsi 


0.007,0.0065 (measured) 


Momentum spread 


A/j 

P 


0.001 (measured) 


Lattice parameter 


Px,tripi Py,trip 


5.49, 7.76 



E. Beam parameters during the measurements 

Experiments were carried out using N 7+ and U 73+ ion beams at the SIS-18 injection energy of 
11.4MeV/u. The data were taken during 600 ms long plateaus. At injection energy space charge 
effects are usually strongest. Four bunches are formed from the initially coasting beam during adi- 
abatic RF capture. The experiment was repeated for different injection currents. At each intensity 
level several measurements were performed with different types and levels of beam excitation in 
both planes. Tune measurements were done simultaneously using the TOPOS and BBQ systems. 
The beam current and the transverse beam profile are measured using the beam current trans- 
former |19j and the ionization profile monitor (IPM) [20] respectively. An examplary transverse 



beam profile is shown in Fig. 13 . The dipole synchrotron tune (Qsi) is deduced using the residual 
longitudinal dipole fluctuations of the bunches. Q s \ has been used as an effective synchrotron tune 
for all experimental results and will be referred as Q s from hereon. The momentum spread is 
obtained from longitudinal Schottky measurements [2T]. Important beam parameters during the 



experiment are given in Tab. I and Tab. II . It is important to note that all the parameters 
required for analytical determination of q sc , q c are recorded during the experiments. 

From |Tab. I and Tab. II one can estimate that in the measurements the head-tail space 



charge and image current parameters were in the range q sc < 10 and q c < 0.2q sc for the horizontal 
and vertical planes. 
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TABLE II. Beam parameters during the N 7+ experiment 



Beam/Machine Parameters 


Symbols 


Values 


Atomic mass 


A 


14 


Charge state 


q 


7 


Kinetic energy 


Ehin 


11.56 MeV/u (measured) 


Number of particles 


N p 


3, 6, 11, 15 • 10 9 (measured) 


Tune 


Qx 7 Qy 


4.16, 3.27 (set value) 


Chromaticity 


£,x i £y 


-0.94, -1.85 (set value) 


Transverse emittance 


e x ,e y (2a) 


33, 12 mm-mrad (measured) 


Slip factor 


V 


0.94 


Bunching factor 


Bf 


0.37 (measured) 


Synchrotron tune 


QsOi Qsi 


0.006,0.0057 (measured) 


Momentum spread 


P 


0.0015 (measured) 


Lattice parameter 


Px,tripi Py,trip 


5.49, 7.76 




FIG. 13. Normalized transverse beam profile in vertical plane at N 7+ = 15 • 10 8 ions. The dotted lines 
shows the normal distribution for the rms width obtained by evaluation of beam profile around its centre. 



V. EXPERIMENTAL OBSERVATIONS 



Tune spectra measurements at different currents are presented and interpreted in comparison 



with the predictions of Sec. Ill . The effect of different excitation types and power on the transverse 
tune spectra is studied. Transverse impedances in both planes are obtained from the coherent tune 
shifts. Incoherent tune shifts are obtained from the relative frequency shift of head-tail modes in 
accordance with |Eq. 8 . Chromaticity is measured using head-tail eigenmodes and the obtained 
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relative height of the observed peaks for different chromaticities is analyzed. 



A. Modification of the tune spectrum with intensity 




FIG. 14. Horizontal tune spectra for t/ 73+ ions and beam parameters given in Tab. I (see text). The 
dashed lines indicate the head-tail tune shifts from |Eq. 8| . 



Figure 14 shows the horizontal tune spectra obtained with the BBQ system using band-width 



limited noise at different intensities. Figure 14 (a) shows the horizontal tune spectrum at low 
intensity. Here the k = 1,0,-1 peaks are almost equidistant, which is expected for low intensity 
bunches. The space charge parameter obtained using the beam parameters and Eq. 4 is q sc ~ 0.15. 



The vertical lines indicate the positions of the synchrotron satellites obtained from Eq. 8| (with 



Figure 14 (b) shows the tune spectrum at moderate intensity (q sc ~ 0.7). The 



k 



2, —2 peaks can both still be identified. Figure 14 (c) shows the tune spectra at larger 



intensity (q sc ~ 1.7). An additional peak appears between the k = and k = — 1 peaks which can 
be attributed to the mixing product of diode detectors (since at this intensity 30 — 40V^ acts across 
the diodes pushing it into the non-linear regime). The k = 0, 1, 2 peaks can be identified very well, 



whereas the amplitudes of the lines for negative k already start to decrease (see Sec. Ill ). In the 
horizontal plane the effect of the pipe impedance and the corresponding coherent tune shift can 
usually be neglected because of the larger pipe diameter. 



Figure 15 shows the vertical tune spectrum obtained by the BBQ system with band limited 



noise excitation for N + beams with q sc values larger than 2. Here the negative modes (k < 0) 
could not be resolved anymore. In the vertical plane the coherent tune shift is larger due to the 
smaller SIS-18 beam pipe diameter (q c & q sc /W). The shift of the k = peak due to the effect of 
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FIG. 15. Vertical tune spectra for N 7+ ions and beam parameters given in Tab. II . The clashed lines 
indicate the head-tail tune shifts from |Eq. 8| . 



the pipe impedance is clearly visible in |Fig. 15 



In the measurements the width of the peaks is determined by the cumulative effect of non-linear 
synchrotron motion, non-linearities of the optical elements, closed orbit distortion, tune fluctuation 



during the measurement interval as well as due to the intrinsic Landau damping ( Sec. Ill ). From 



the comparison to the simulations we conclude that the intrinsic Landau damping is an important 
contribution to the width of the k = 1, 2 peaks. 

B. Determination of coherent and incoherent tune shifts 



The coherent tune shift AQ C can be obtained by measuring shift of the k = line as a function 



of the peak bunch current as shown in Fig. 16 . The transverse impedance is obtained by a linear 



least square error fit of the measured shifts in both planes to Eq. 6 and Eq. 7 . The impedance 
values are obtained in the horizontal and vertical planes at injection energy are 0.23 MO/m 2 and 
1.8 Mfl/m 2 respectively, which agrees very well with the expected values for the average beam 
pipe radii of the SIS- 18. 



Figures 17 18, 19 and 20 show the measured positions of the peaks in the tune spectra for 



different intensities. In comparison the analytical curves (solid lines) obtained from Eq. 8 for 



the head-tail tune shifts are plotted using q sc estimated from the beam parameters in [Tab . I 
for each intensity. The error bars in the vertical plane (5qk = ^p 1 ) correspond to 



and 



Tab. II 



the 3 dB width of the measured peak due to accumulation of various effects (see Sec. VA ). In 



the horizontal plane, error bars (5q sc ) are estimated by propagation of parameter uncertainties 
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FIG. 16. Coherent tune shift obtained from the measurement for the horizontal and for the vertical planes 
as a function of the peak beam current. The dotted lines correspond to a linear least square error fit. The 
error bars in the horizontal plane are due to uncertainties in the current measurements. In the vertical plane 
errors result from the width of k = mode. The FFT resolution is 5 • 10 -4 and is always kept higher 
than the mode width. 




FIG. 17. The measured positions of the peaks in the horizontal tune spectra for different U 73+ beam 



intensities together with the analytical curves from Eq. 8 using the space charge tune shift estimated from 



the beam parameters in Tab. I . The dotted line corresponds to the incoherent tunc shift. The error bars 
for the vertical plane correspond to the width of measured modes (see text). In the horizontal plane the 
error bars are estimated by the propagation of uncertainties (see text). 



mentioned in ISec. VIll . 

In this subsection, we introduce another space charge parameter q sc ,m which is the measured 
space charge parameter using the following method. It is not to be confused with q sc which is 



predicted for a given set of beam parameters by Eq. 4 . The incoherent space charge tune shift 
can be determined directly from the tune spectra by measuring the separation between the k = 
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FIG. 18. This plot shows the predicted shifts from analytical |Eq. 8| and measured head-tail mode frequencies 



are overlaid in vertical plane for U 73+ ion beam at various current levels ( Tab. I ). 




FIG. 19. This plot shows the predicted shifts from analytical |Eq. 8| and measured head-tail mode frequencies 
are overlaid in horizontal plane for N 7+ ion beam at various current levels ( 



Tab. II 



and k = 1 peaks, i.e. (qkoi = A g fc01 ) and fitting it with the parameter q sc in the predictions 



from Eq. 8 . The value of q sc for the best fit is denoted as q s 



1 



fffcOl 

7%oi 



< %oi < 1 
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Equation 19 is obtained by rearranging Eq. 8 for k 
absolute error on measured q scm (Sq S c,m) is given by 



0,1 while 7 = qsc „ 9c . The linearized 

Qsc 



Sq sc ,r 



-(1 + gfcOl) 



IkOl 
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FIG. 20. This plot shows the predicted shifts from analytical |Eq. 8| and measured head-tail mode frequencies 



are overlaid in vertical plane for N 7+ ion beam at various current levels ( Tab. II ) 



fiQkoi is given by either the width of the k = 0, 1 lines or by the frequency resolution of the system. 
In a typical tune spectra measurement using data from 4000 turns dqkoi ~ 0.04. The absolute 



error is a non- linear function of q sc in accordance to the Eq. 19 . It is possible to define the upper 
limit of q sc where this method is still adequate based on the system resolution and |Eq. 19 . If 
we define a criterion that, qkoi > Squoi to resolve the head-tail modes. This gives the limit to be 
Qsc 8 where the measurement error is still within the defined criterion. 



Figure 21 shows a plot of the predicted space charge tune shifts (q S c) versus the ones measured 



from the tune spectra using the above procedure (q sc ,m)- For q sc < 3.5 the space charge tune shifts 
measured from the tune spectra are systematically lower by a factor = 0.74 than the predicted 



shifts. It is shown by the dotted line in Fig. 21| which is obtained by total least squares fit of the 
measured data points. For larger q sc the factor decreases to ~ 0.4. Thus the method for measuring 
the incoherent tune shift based on head-tail tune shifts is found to be satisfactory only in the range 
qsc ;$ 3.5. A possible explanation is the effect of the pipe impedance. Similar observations are 



made by the results of self-consistent simulations in Sec. Ill , where for q sc > 2 the separation of 



the k = and k = 1 peaks observed is underestimated by Eq. 8 (see Fig. 6 ). 



C. Effect of excitation parameters on tune spectrum 



Figure 22 presents the tune spectra obtained from BBQ system at various excitation power 
levels of band limited noise. The beam is excited with 0.25, 1.0 and 2.25 mW/Hz power spectral 
density on a bandwidth of 10 KHz. Signal-to-noise ratio (SNR) increases with excitation power 
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FIG. 21. Combining the results from the Figures [TTJ [18] [19] and [20} a plot of predicted q sc using |Eq. 4| 
against measured g sc using the distance between modes k = and fc = 1 is obtained. 

whereas the spectral position of various modes is independent of excitation power. 
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FIG. 22. Tune spectra at 5 • 10 8 f7 73+ particles with beam excited at three different excitation amplitudes. 
The beam is excited with 0.25, 1.0 and 2.25 mW/Hz power (blue, red and black respectively). The dotted 
lines mark the relative positions of the head-tail modes. The frequencies of various head-tail modes are 
unaffected. 

Beam excitation using two other excitation types i.e. frequency sweep and white noise is also 



performed to study the effect of excitation type on the tune spectra. Figure 23 shows the tune 
spectra under same beam conditions for different types of beam excitation obtained from the 
BBQ system. The frequencies of various modes in the tune spectra are independent of the type of 
excitation. The signal-to-noise ratio (SNR) is optimum for band limited noise due to long averaging 
time compared to "one shot" spectra from sweep excitation. 
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D. Time domain identification of head tail modes 



Figure 24 shows the 2-D contour plot for frequency sweep excitation in vertical plane obtained 



from the TOPOS system, where various head-tail modes are individually excited as the excitation 
frequency crosses them. Frequency sweep excitation allows resolving the transverse center of mass 
along the bunch for various modes which helps in identifying each head-tail mode in time domain. 
This serves as a direct cross-check for the spectral information and leaves no ambiguity in identi- 



fication of the order (k) of the modes. Figure 25 shows the corresponding transverse center of 
mass along the bunch for k = 0, 1 and 2 at the excited time instances. This method works only 
with sweep excitation and requires high signal-to-noise ratio in the time domain, which amounts 
to higher beam current or high excitation power. 




FIG. 23. Horizontal tune spectra over time with U 73+ at 5 TO 8 particles with band limited noise (RF Noise), 
frequency sweep (chirp noise) and white noise (top to bottom). The dotted lines mark the relative positions 
of the head-tail modes. 



E. Measurement of the chromaticity 



As highlighted in the previous section, the frequency sweep allows to resolve the different head- 
tail modes both spectrally and temporally. This procedure can be used for the precise determination 



of the chromaticity by fitting the analytical expression for the head-tail eigenfunction Eq. 9 to 



the measured bunch offset, with the chromaticity (£) as the fit parameter as shown in Fig. 26 for 
k = 0, 1 and 2. The measured chromaticity is independent of the order of head-tail eigenfunction 
used to estimate it. 
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FIG. 24. Vertical tune spectra with frequency sweep over time at N 7+ with 1.5 • 10 9 particles. The various 
modes get individually excited as the sweep frequency (marked by the arrow) coincides with the mode 
frequency. The k — 0, 1, 2 modes arc marked. Symmetric sidebands due to coherent synchrotron oscillations 
are visible around the sweep frequency spaced with Q s and are clearly distinguished from the head-tail 
modes. 
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FIG. 25. Transverse center of mass along the bunch for k = 0, 1 and 2 modes corresponding to |Fig. 24| . 
Frequency sweep is used to excite the beam which enables to resolve the distinct head-tail modes temporally. 
The rms bunch length (2ay ms ) is indicated. 



The fitting method is shown in |Eq. 21 ; the head-tail eigenfunction from Eq. 9 is multiplied 



with the beam charge profile S(t) and corrected for the beam offset Ax at the BPM where the 
signal is measured. 



F(4>, £, A) = Ax ■ S(t) ■ (1 + A - x{4>)) 
£(^) = (M(0)-F(<M)) 2 
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The fit error E(£, A) is reduced as a function of independent variables; chromaticity £ and head-tail 
mode amplitude A. The fit error gives the goodness of the fit. It is used to determine the error bars 
on measured chromaticity. This method has been utilized for the determination of chromaticity at 



SIS-18 as shown in Fig. 27 . The set and the measured chromaticity can be fitted by linear least 
squares to obtain the form £ S)3/ = 1.187£ TO>2/ + 0.804 as shown by red dashed line in Fig. 27 . Fig- 



ure 27 also shows a coherent tune shift due to change in sextupole strength which is used to adjust 



the chromaticity. This is due to uncorrected orbit distortions during these measurements. These 
chromaticity measurements agree with the previous chromaticity measurements using conventional 
methods 



It is also possible to determine the relative response amplitude of each head-tail mode to the 



beam excitation both in time and in frequency domain with TOPOS. Figure 28 shows the tune 
spectrum obtained with sweep excitation for different chromaticity values. The beam parameters 
are kept the same (N 7+ , 14 • 10 9 ,q sc ~ 10). The spectral position and and relative amplitude of 



each head-tail mode peak are confirmed using the time domain information (see Fig. 25 ). In 



Fig. 29 the single particle response amplitudes for different k ( Eq. 3 ) are plotted as a function 



of the chromaticity. The measured relative amplitudes are indicated by the colored symbols. 



The comparison indicates that the simple single particle result ( Eq. 3| ) describes quite well the 
dependence of the relative height of the peaks obtained from the TOPOS measurement. 




FIG. 26. Analytical curves from Eq. 9 (dotted) are fitted to the local transverse offset for the k = and 
k = 1 modes with the head-tail phase shift as the fit parameter. This method has been used for the precise 
determination of the chromaticity. 



25 



O 0.3 
% 0.28 

CO 

2 0.26 



-3.5 



-2.5 -2 -1.5 
Measured E 



-0.5 



FIG. 27. Chromaticity measurement using the method shown in Fig. 25 over the possible range of operation 
at SIS-18. The black dashed line marks the measured natural vertical chromaticity of SIS-18. The error 
bars for the chromaticity measurement are determined mainly by the electronic noise in the amplifier chain 
of the TOPOS system. 
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FIG. 28. Relative signal amplitudes of various head-tail modes are plotted with respect to the measured 
chromaticity. The centroid of the spectrum moves toward higher order modes with increase in chromaticity. 



VI. APPLICATION TO TUNE MEASUREMENTS IN SIS-18 



In this section we will discuss the application of our results to tune measurements in the SIS-18 
and in the projected SIS-100, as part of the FAIR project at GSI [21]. As shown in the previous 
section the relative amplitudes of the synchrotron satellites in the tune spectra are primarily a 
function of chromaticity and possibly the excitation mechanism. In order to determine the coherent 
tune with high precision the position of the k = mode has to be measured. Depending on the 
machine settings, if the relative height of the k = peak with respect to the other modes is small, 
then the k = mode may not be visible at all. To estimate the bare tune frequency in this case, the 
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FIG. 29. Relative signal amplitudes of various head-tail modes are plotted with respect to the coefficient of 



Bessel's function in Eq. 3 (xf = x/2)- All the modes are relative to k=I mode, which has a significant 
amplitude during the whole measurement range. 



information of space charge parameter, coherent tune shift and chromaticity are all simultaneously 
required with good precision. 

Another important point is the tune measurement during acceleration. The space charge pa- 
rameter for 1 • 10 10 Ar 18+ stored ions in the SIS-18 from injection to extraction reduces only by 



~ 20% as shown in Fig. 30 . The dynamic shift of head-tail modes during acceleration is shown 
in 



Fig. 31 obtained from TOPOS system under same conditions. The asymmetry of k = 1,-1 



modes around k = mode can only be understood in view of the space charge effects predicted 



by Eq. 8 . Thus, a correct estimate of this parameter plays an important role in understanding the 
tune spectra not only during dedicated experiments on injection plateau, but also during regular 
operations. 

The measurement time required to resolve the various head-tail modes (AQ&) is a complex 
function of Q s , q S c, beam intensity and excitation power. To give some typical numbers for SIS-18; 
on a measurement time of 600 ms on the injection plateau, if one spectrum is obtained in ~ 20 
ms (~ 4000 turns), an improvement of factor ~ 6 in SNR by averaging 30 spectra. Following the 



calculations in Sec. VB , q sc < 8 can be resolved under typical injection operations. However, the 
constraints on measurement time are much higher during acceleration, where the tune/revolution 
frequency increases due to acceleration. This allows the measurement of a single spectrum typically 
only over 500 — 1000 turns (depends on ramp rate as well). There are no averaging possibilities 
since the tune is moving during acceleration due to dynamic changes in machine settings as seen 
in Fig. 31 . In addition, the synchrotron tune reduces with acceleration making it practically very 
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difficult to resolve the fine structure of the head-tail modes for q sc > 2. 
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FIG. 30. Change in synchrotron tune and space charge tune shift on acceleration from injection to extraction 
for Ar ls+ with 1 • 10 10 particles. The estimated space charge parameter reduces by w 20% in this typical 
case from injection at 11.4 MeV/u to 300 MeV/u at extraction. 
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FIG. 31. The movement of head-tail modes during acceleration. The head-tail modes are getting closer 
since the synchrotron tune reduces with acceleration, but the asymmetry of the k — 1 and k = — 1 modes 
around k = is maintained throughout the ramp which depends primarily on the space charge parameter 



VII. CONCLUSION 



Two complimentary tune measurement systems, TOPOS and BBQ, installed in the SIS-18 are 
presented. Analytical as well as simulation models predict a characteristic modification of the tune 
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spectra due to space charge and image current effects in intense bunches. The position of the 
synchrotron satellites corresponds to the head-tail tune shifts and depends on the incoherent and 
the coherent tune shifts. The modification of the tune spectra for different bunch intensities has 
been observed in the SIS- 18 at injection energy, using the TOPOS and BBQ systems. From the 
measured spectra the coherent and incoherent tune shift for bunched beams in SIS- 18 at injection 
energies were obtained experimentally using the analytic expression for the head-tail tune shifts. 
Head-tail modes were individually excited and identified in time domain and correlated with the 
spectral information. A novel method for determination of chromaticity based on gated excitation 
of individual head tail mode is shown. The dependence of the relative amplitudes of various head- 
tail modes on chromaticity is also studied. The systems were compared against each other as well 
as with different kinds of excitation mechanisms and their respective powers. These measurements 
give a clear interpretation of tune spectra at all stages during acceleration under typical operating 
conditions. The understanding to tune spectra provides an important input to new developments 
related to planned transverse feedback systems for SIS-18 and SIS-100. The measurement systems 
also open new possibilities for detailed beam investigations as demonstrated in this contribution. 
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APPENDIX 



Calculation of bunching factor 



Figure 32 shows a typical longitudinal profile of the bunch. If Sj is the amplitude at each time 



instance j = 1, .., iV. The bunching factor is calculated by the Eq. Al 



B 



f 



Al 



max(Sj) ■ N 

where N ~ 147 is the number of samples in one RF period (at injection). The TOPOS system 
samples the bunch at 125 MSa/s, thus the difference between adjacent samples is 8ns. 
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FIG. 32. A typical normalized longitudinal beam profile from TOPOS system at injection. The length of 
the RF period is w 1.2 [is. 

Measurement uncertainties 



The calculation of AQ SC from |Eq. 4 has a dependence on measured current, transverse beam 
profiles, longitudinal beam profiles and the twiss parameters. The measurement uncertainty on 
each of these measured parameters at GSI SIS- 18 were commented in the detailed analysis in [23J. 
Even though some parameters and the associated uncertainties are correlated, any correlations are 
neglected in present analysis. Uncertainties in each measured parameter are propagated to find 
the error bars on the calculated and measured incoherent tune shifts. 

Reproducing from Ref. [23], the relative random uncertainty (std. deviation) in beam profile 



width(o" x ) measurements is given by Eq. A2 



5a x 0.043 



+ 0.335 Ni A2 



a x B 

where Ax = 2.1mm is the wire spacing and 5N{ = ^ is the ADC resolution of the IPM and 
B is defined as a x / Ax. If the error bars are derived from j measurements, the measured profile is 
given by 

cr a v,x = {&x)j, Sa av ,x = J {ofjj - {<Tx) 2 j + (S(Tx) 2 A3 
For each tune measurement at the given intensity and excitation power, 5-8 transverse beam profiles 
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were measured, and the relative error is obtained ~ 5% using Eq. A2 and Eq. A3 . The relative 
systematic error (bias) in transverse beam width measurements is < 1% [23] and ignored in this 
analysis. 

The uncertainty in the injected current is dominated by fluctuations in the source and the 
relative uncertainty is estimated to be ~ 5% based on 5-8 measurements at the same intensity 
settings for each measurement point. Bunch length and bunching factor vary by ~ 2 — 3% due to 
long term beam losses only under high intensity beam conditions. The maximum relative bias in 
the lattice parameter B is assumed to be ~ 5% at the IPM location. Taking all the relative errors, 



uncertainty propagation using familiar Eq. A4 gives relative error for estimated incoherent tune 
shifts pa 12%. 
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Tune measurements done by averaging over long intervals contribute to the width of modes due 
to long term beam losses. Beam losses lead to change in coherent tune especially in the vertical 
plane where the image current effects are larger. This has been highlighted at appropriate sections 
in the text. 
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